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Abstract 

Polyakov has suggested that two dimensional turbulence might be described by a min- 
imal model of conformal field theory. However, there are many minimal models satisfying 
the same physical inputs as Polyakov's solution .M(2,2i)- Dynamical magnetic fields and 
passive scalars pose different physical requirements. For large magnetic Reynolds number 
other minimal models arise. The simplest one, A^(2,i3), makes reasonable predictions that 
may be tested in the astrophysical context. In particular, the equipartition theorem be- 
tween magnetic and kinetic energies does not hold: the magnetic one dominates at larger 
distances. 
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Turbulence is one of the longstanding problems of theoretical physics. It affects an 
enormous variety of phenomena, from galactic plasma to cigarette smoke [1]. As such, 
the problem is far from being solved. The randomness of the basic physical quantities 
clearly suggests that the theory of random fields should be used to study this problem. 
However, since one is dealing with a strongly coupled field theory problem, the lack of 
reliable computational techniques has slowed down the progress. 

Recently, Polyakov [2] has made a bold attempt to understand fully developed, homo- 
geneous, isotropic, stationary two dimensional turbulence by using conformal field theory. 
Fully developed turbulence simply means that the motion of the fluid is not affected by ex- 
ternal constraints such as boundaries (at large distances), or viscosity (at short distances). 
It is in this case that universal features of turbulence have a better chance to be under- 
stood. For fully developed turbulence, the assumptions of spatial homogeneity (translation 
invariance) and isotropy (rotation invariance) are very reasonable. More subtle is the sta- 
tionary assumption (time translation invariance) because it implies that turbulence must 
be sustained by external forces. The restriction to the two dimensional case is, of course, 
because one wants to use the full power of the conformal field theory [3] . 

Polyakov's basic idea is to describe the universal features of two dimensional turbu- 
lence by fitting it into a (non unitary) minimal model M( p , q ) of two dimensional conformal 
field theory, here 1 < p < q, p and q integers, prime to each other [3]. Any such model 
contains (p-l)(q-l)/2 primary fields $( m , n ), < m < p, < n < q, $( m , n) = $( p _ mi9 _ n ) 
with conformal weight 

(pn - qm) 2 - (p - q) 2 

\m,n) = ^ • (1) 

For completeness, let us recall that the fusion rules are 

min ( m ^ -\-m>2 — 1 , rain(n ^ +^2 — ^ > 
2p — m.-^ — m.2 — 1) 2q — n-^ — tt-2 - 1) 

^(rm.m)] >< [$(m 2l » 2 )] = J2 ' J2' ^,1)], (2) 

k=\mi— m 2 | + l l=\ni— n 2 | + l 

where the symbol £' means that the sum is performed with index jumping of two by two, 
e.g. k = 3,5,7 - • •. The infrared contributions from the regular terms must be always 
kept in mind because otherwise exact conformal symmetry would imply the orthogonality 
theorem between primary fields. 



In two dimensions, the velocity field v a of an incompressible fluid {d a v a = 0), can be 
written in terms of a single (pseudo) scalar function ip, called the stream function: 

v a = eapdpip. (3) 

The vorticity uj associated to the velocity field is a scalar and can also be written purely 
in terms on the stream function: 

w = £apd a vp = -d a d a ip. (4) 

One then attempts to interpret tp as one of the primary fields of some minimal model 
ip = $( m , n ) and to use the dynamics of the theory (Navier-Stokes equations plus energy 
consideration) to fix p, q, m and n. Let us briefly recall how this is done in [2]. 

The Navier-Stokes equations for this two dimensional incompressible fluid can be 
written, (after taking the curl) as 

= dj + e a pd^d a uj = ud a d a ip. (5) 

The notation is as follows: throughout this paper, the symbol ^ will denote the convective 
derivative of any quantity, i.e., the time derivative along the motion of a particle in the 
fluid. The ordinary partial time derivative will be denoted by a dot or by d t . The relation 
between the two is 

= d t + v a d a = d t + e^dpipda. (6) 

The constant v represents the kinematic viscosity. In the following we will always work in 
the so called inertial range, where the effective Reynolds number is huge and the viscous 
term in (5) can be neglected; 

I" = °- (7) 

One then talks of (7) as the inviscid Navier-Stokes equations. 

In two dimensions, in addition to the kinetic energy per unit mass density E = 
\ f v a v a d 2 x, there is another quadratic conserved quantity. This new quantity is defined 
as the integral of the square vorticity per unit mass density: H = | J u> 2 d 2 x and it is 
known as the "enstrophy" . From the detailed analysis in [4] it turns out that, for scales 



smaller than the scale at which energy is fed into the system, the flux of enstrophy rj is 
the only relevant dimensional parameter and one could estimate the scaling behavior of all 
observables in the Kolmogorov approximation. For example, the energy spectrum E(k), 
defined as 

^ < VaVa >= J d 1 kE(k), k = y/k a k a , (8) 

is given, by dimensional analysis, as E(k) ~ r] 2 ^ s k~ 3 . Of course, the aim is to find 
corrections to this scaling behavior. From the definition of E{k) and being v a a level one 
spin one field in [ip] one has 

E(k) » A;- 1+2(A « +A ^ = k 1+4A ^ (9) 

and what is left is to find A^. 

In order to achieve this goal, one first translates the constant flux of enstrophy into 
the scaling behavior of the two point correlation function: 

< dj(r)u(0) >^ r°. (10) 

Now, if the primary field ip satisfies the fusion rules 

[ip] x [ip] = [(f)] + less singular terms, (11) 

and 4> is neither the identity nor a level two degenerate field, from (7) one has uj E [4>] at 
level two and from (4) one has u G [ip] at level one. Hence: 

= A w = + 1 and A^ = A^ = A^ + 2. (12) 

Equation (10) then can be used to fix the sum of the conformal weights: 

A^ + A^ + A„ + A w = 2(A^ + 2 + A^, + l) = 0. (13) 

The physical conditions coming from the inviscid Navier-Stokes equations and the 
enstrophy flux then prompt us to look for a minimal model M.( p , q ) containing two primary 
fields satisfying the following two conditions: 

[ip] x [ip] = [</>] + less singular terms; and A<^ + = —3. (14) 
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Polyakov has found the following solution: M^,2i)i with ip = $(1,4) and <fi = $(1,7)- Un- 
fortunately, there are many (perhaps infinitely many) other solutions to these conditions. 
A very detailed list of solutions with energy spectrum steeper than k~ 3 is given very re- 
cently in [5] . We also independently found the first few of them with the help of a simple 
computer program. In the table below we summarize the first few solutions we have found, 
with no particular constraints on the energy spectrum. The first one is Polyakov's solu- 
tion. Together with his solution, ^(3,25), M-(3,26) an< 3 -M-(6,55) gi ye the most reasonable 
energy spectrum, i.e., steeper than k~ s . These are among those found by Matsuo in [5], 
where the condition on the steepness of the spectrum is enforced by the extra constraint 
< —1. Notice that only the first and the last solutions are parity invariant, in the 
sense [ip] £ [ip] x [ip] , but there are others, not included in the table. 

First few solutions of Polyakov's requirements 



(p, q) 




<t> 


E(k) « 


(2,21) 


$(1,4) 


$(1,7) 


£.-25/7 


(3,25) 


$(1,11) 


$(1,9) 


£-23/5 


(3,26) 


$(1,5) 


$(1,9) 


£-55/13 


(3,35) 


$(1,21) 


$(1,11) 


k~ 9 / 7 


(5,51) 


$(1,17) 


$(1,11) 


£-47/17 


(6,55) 


$(1,14) 


$(1,9) 


£-41/11 



These results show that the requirements above are not restrictive enough to pin 
down uniquely one minimal model. Nevertheless, Polyakov's solution is somewhat special 
in that it contains the smallest number of primary fields and hence might be more stable. 
It is interesting to include magnetic fields into the picture to see if constraints are more 
stringent. But before we discuss this attempt, let us first consider a simpler situation: the 
diffusion of a passive scalar into Polyakov's M(2,2i) solution. 

The term "passive scalar" denotes any scalar field such as temperature, density of a 
pollutant etc., that does not modify the Navier-Stokes equations and satisfies the diffusion 

5 



equation 

jT = Kd a d a T. (15) 

The assumption that T does not appear into the Navier-Stokes equations always implies 
some approximation. For instance, if T stands for the absolute temperature, the assump- 
tion is that the dilation coefficient of the fluid is small, i.e., the density remains constant 
over the temperature range being examined. From now on we will identify T with the 
absolute temperature. 

We are interested in the so called inertial-convective range, where the Navier-Stokes 
equation is (7) and the equation for the temperature reduces to = 0. The temperature 
spectrum is defined, in analogy with the energy spectrum, as 

^ < T 2 >= J dP-kErik). (16) 

Let us first look at the spectrum in the Kolmogorov approximation. Let us define the 
temperature dissipation rate as 6t- Clearly, 

dim(E T ) = dim(T 2 ) x length and dim(e T ) = dim(T 2 ) x time -1 (17) 

The temperature spectrum must then be proportional to €t and, since in the enstrophy 
dominated inertial range the only other dimensional parameter is rj, one obtains the well 
known scaling law 

Erik) « eTV^^k' 1 (18) 

This is the Kolmogorov-type prediction of the spectrum in the inertial-convective range. 
Note, in particular, that it is obeyed by the vorticity u itself when treated as a passive 
scalar, being E UJ (k) = k 2 E(k) and e u = rj. 

The universality of this behavior suggests that, in this regime, the exponent for the 
temperature spectrum should be the same as that for the vorticity spectrum, even if the 
latter is not simply given by the Kolmogorov exponent —1 above. The M^,2i) solution 
predicts Ex(k) ~ /c -11 / 7 . In any case, it is not correct to try to identify another primary 
field in as the temperature because the predicted exponents would be very far off. 

Let us now turn to the main focus of this paper: two dimensional turbulence in pres- 
ence of magnetic fields. In three dimensional space, the magnetohydro dynamics (MHD) 
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equations can be written, in the limit of infinitely conducting, incompressible and inviscid 

— * 

fluid, purely as functions of the velocity v and the magnetic field B. It is convenient to 
rescale the magnetic field as given in S.I. units by the constant factor l/yfjlp, so that the 
rescaled field has the dimensions of velocity. With this convention, taking the curl of the 
Navier-Stokes equation to get rid of the pressure term, the equations take the form 

V -v = 
V-B = 

V - [ 0,'r) = V x ( - (r- V)r+ (B ■ V)B ' 



d t B = (B-V)v- {v-V)B 

Let us now consider the case of an "effectively two dimensional" theory, defined by 
the condition: 

d 3 v = and d 3 B = 0. (20) 

With this assumption, one can treat the third components of the magnetic field and of the 
velocity as a two dimensional scalar: B 3 = B and = v. Furthermore, one can introduce 

— * _ 

two stream functions A and ip for the first two components of B and v respectively: 

B a = tapdpA and v a = e a pdpij). (21) 

For this effectively two dimensional theory, equations (19) can be written simply in terms of 
four independent two dimensional scalars: A, B, v and ip. These equations look particularly 
nice if one also introduces (as dependent quantities) the two dimensional curls 

J = tapdaBp = -d a d a A and u> = e a pd a vp = -d a d a ip, (22) 

and the differential operator: 

A = e aP dpAd a . (23) 
With the aid of the quantities defined above, the equations of motion can than be written 
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as 



d 



AJ 




AB 



Av 



(24) 



0. 



This shows that the fields B and v only appear in the second and the third equations of 
(24). It is therefore consistent to study first the possible turbulent solution of A and by 
setting B and v to zero. 

Before we go any further, one important remark must be made, in order to avoid 
confusion. The kind of turbulent solution we are considering is not the same as the kind of 
two dimensional turbulence that has recently been studied in laboratory [6] . In laboratory 
experiments, one is dealing with magnetic Reynolds numbers that are fairly small and the 
system (typically a tank filled with mercury) is subjected to a strong external magnetic 
field perpendicular to the plane. One then writes the total magnetic field as the sum of 

— * — * 

the external one -B ex t. and the turbulent one b. Since the turbulent component is much 
smaller than the external field, the MHD equations can be expanded in first order in b. 
After taking the curl, the magnetic field drops out completely of the Navier-Stokes equation 
and its only effect is to reduce the B = three dimensional Navier-Stokes equation to the 
B = two dimensional Navier-Stokes equation (5). One could then study the behavior 

— * 

of b as a passive vector on the ordinary turbulent background. As far as conformal field 
theory goes, one is then brought back to Polyakov's solution. 

In the following, we consider a different regime altogether. We assume that there is 
no external magnetic field and that the turbulent magnetic field is parallel to the plane 
and is (in our units) comparable with the velocity. This can happen only at very high 
magnetic Reynolds number, such as in astrophysical situations (for instance in a galaxy). 
One last word of caution: we are not interested in the dynamo effect. As a matter of fact, 
the dynamo effect does not take place at all in two dimensions [7]. In two dimensions, 
in order to keep turbulence in the stationary regime, one has to pump in magnetic energy 
as well as kinetic energy, otherwise the magnetic field would dissipate due to the finite 
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conductivity of the material just as much as the velocity would dissipate due to the finite 
viscosity. 

The third component of the magnetic field need not be ignored but can be treated 
at a later stage since it decouples from the equations for ip and A as already stressed. In 
any case, one cannot invoke this magnetic field as the agent that blocks the system into 
the two dimensional state. However, if one is willing to take the astrophysical scenario 
seriously, there are other possible effects that might help, such as the rotation of the galaxy 
or gravity itself. 

Let us start by rewriting the equations for ip and A more explicitly 

w + e a pdf3il)d a Lo = e a pdpAd a J 

(25) 

A + eapdpipdaA = 0. 

From (25) we see that a radical change takes place in this regime. Enstrophy is no longer a 
conserved quantity. Its place is taken by the integral of the square of the stream function 
A. The only two quadratic conserved quantities are then 

E =\J ( v *v a + B a B a )d 2 x and G=\J A 2 d 2 x. (26) 

It is then natural to look for a solution for which 

< A(r)A(0) >^ r°. (27) 

Instead of the constant enstrophy condition, we are considering €a, the dissipation rate of 
A, as a constant. 

Let us translate these physical conditions in the language of conformal field theory. 
Again, the working assumption is that we try to fit the requirements (25) and (27) by a 
minimal model. We interpret the two basic scalars A and ip as primary fields and A as 
a secondary field in the conformal family [x] of the most singular field appearing in the 
fusion rules of A and ip. By using point split regularization, one can see that A must be a 
level two scalar in [%]. 
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We must then look for a minimal model Ai( P;q ) satisfying the following requirements: 

a) There are two primary fields A and x f° r which + A x = —2. 

b) There is another primary field ip such that its fusion rules are, up to less singular 
terms, [ip] x [A] = [%]. 

To be sure, these conditions do not fix the theory uniquely. The simplest solution 
however, provides a nice example. Surprisingly, this solution is simpler (i.e., contains less 
primary fields) than Polyakov's. It is the minimal model -M(2,i3), whose basic structure is 
summarized in the table below. 

Structure and interpretation of -M( 2 ,i3) 





$(i,i) 


$(1,2) 


$ (1,3) 


$(1,4) 


$(1,5) 


$ (1,6) 


meaning: 


I 




? 


A 


X 


? 


A = 





-5/13 


-9/13 


-12/13 


-14/13 


-15/13 



It should not be surprising that this solution is not the same as the one before. Remember 
that the introduction of magnetic fields adds an extra dimensional parameter to the theory. 
This parameter is still relevant in the inertial range and it affects the scaling laws already 
in Kolmogorov theory. 

Having found a solution for our physical requirements, it is now a simple exercise in 
conformal field theory to work out its basic predictions. The kinetic energy and magnetic 
energy spectra behave respectively as 

E kin .(k) « k~ 7/13 and E mag .(k) ^ k~ 35/ls . (28) 

This means that the equipartition theorem does not hold exactly for our solution. In fact, 
the kinetic energy will dominate at short distance and the magnetic one will take over at 
large distance. 

Passive scalars will still be represented by level one fields. But now, the argument 
given above using the vorticity no longer apply, because u no longer satisfies the equation 
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■^uo = 0. For example, if one takes T e [<E>( 1)6 )], level one scalar, one finds E T (k) ps /c -21 / 13 
but we do not see any physical reason for that identification. 

In conclusion, Polyakov's idea gives a fresh new look on one of the oldest unsolved 
problems of physics. Unfortunately, the physical requirements found so far are not enough 
to uniquely identify one specific minimal model. In fact, there are perhaps infinitely many 
minimal models satisfying the constraints. One wonders whether the "correct" conformal 
field theory might be the one obtained as some limiting case of the minimal models above 
or perhaps as some other, non minimal, conformal field theory. Quite apart from this, the 
introduction of magnetic fields and passive scalars is of some interest. It is still not enough 
to fix the theory completely but the simplest model satisfying the physical requirements 
gives sensible predictions that might be relevant in the astrophysical context. 

We would like to tank S.G. Rajeev for his interest and his helpful comments. We also 
benefitted from discussions with A. Das and P. Bedaque. This work is supported in part 
by DOE Grant DE-FG02-91ER40685. 
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